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The recently suggested quasi-local spin-angular momentum expressions, based on the Bramson superpotential 
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and on the holomorphic or anti-holomorphic spinor fields, are calculated for large spheres near the future null 
infinity of asymptotically flat Einstein-Maxwell spacetimes. It is shown that although the expression based 



, on the anti-holomorphic spinors is finite and unambiguously defined only in the center-of-mass frame (i.e. it 

diverges in general), the corresponding Pauli-Lubanski spin is always finite, free of gauge ambiguities, and 
is built only from the gravitational data. Thus it defines a gravitational spin expression at the future null 
infinity. The construction based on the holomorphic spinors diverges in presence of outgoing gravitational 
radiation. For stationary spacetimes both constructions reduce to the 'standard' expression. 

> 

O ■ !• Introduction 

o> 

The energy-momentum and angular momentum density of matter fields is described by their energy-mom- 
entum tensor T ab , whose contraction with a Killing field K a gives a divergence free current T ab Ki,. The flux 
of this current on a compact, spacelike hypersurface £ with boundary S := is therefore conserved, i.e. 
^ ' depends only on the boundary S and independent of the rest of the hypersurface. In Minkowski spacetime the 

independent Killing 1-forms are the translations K* :— V e x a and the rotations K ah := x a V e x b — x b \7 e x a , 



where x a are the standard Cartesian coordinates, a, b = 0, 3. The corresponding conserved quantities are 
PJ := / s T ab t a K%&Y< and J| b := / E T ab t a K ab dT,, where t a is the future directed unit timelike normal to 



£ and d£ is the induced volume element on S. Pg and J| b can be interpreted as the quasi-local energy- 
momentum and angular momentum of the matter fields associated with the 2-surface S, respectively. If £ 
extends to the spacelike infinity i° or to the future (or past) null infinity I of the Minkowski spacetime, then, 
provided the limits exists (by imposing certain fall-off and global integral conditions, e.g. parity conditions 
at the spacelike infinity), then they define the global energy-momentum and angular momentum. 

In general relativity the usual fall-off conditions on the 3-metric and extrinsic curvature on a spacclikc 
hypersurface ensure that the ADM energy and spatial momentum are finite and unambiguously defined 
[1,2]. However, to have finite and unambiguously defined angular momentum and center-of-mass additional 
conditions, e.g. explicit parity conditions on the three-metric h a b and extrinsic curvature Xab, have to be 
imposed [2,3]. These are preserved by the evolution equations, ensure the functional differentiability of the 
Hamiltonian on the whole phase space, and yield unique Poincare structure both for the lapse-shift at infinity 
and the Hamiltonian. At spatial infinity the energy-momentum and relativistic angular momentum are the 
value of the functionally diffcrcntiablc Hamiltonian, which are just the familiar ADM energy-momentum and 
spatial angular momentum, and the center-of-mass of Beig and O Murchadha [3]. Thus, although there are 
interesting open issues, e.g. whether the angular momentum measured at the spacelike infinity enters the 
Penrose inequality or not [4], the spatial infinity case is well understood. 
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At null infinity there is a generally accepted definition of the energy-momentum, which is the Bondi 
Sachs energy-momentum. However, on the definition of angular momentum there is no consensus at all, 
and there are various suggestions for that: The constructions based on the Komar-Winicour-Tamburino 
expression are intended to be associated with any BMS vector field, and, in the special case of the boost- 
rotation vector fields, they can be interpreted as angular momentum [5-7]. The suggestion of Ashtekar and 
Strcubel is based on symplectic techniques [8], which turned out later to be connected with the Winicour- 
Tamburino linkages [9] . The general form of other expressions are based on charge integrals of the curvature 
[10-16]: In particular, although the suggestion of Penrose [12] is based on the solutions of the 2-surface 
twistor equations and the concept of 'origin' is largely decoupled from the cuts of T + , it fits nicely to the 
symplectic structure of I + [17-20]. Other remarkable suggestion in this class is due to Moreschi [13-16]. He 
considers (and, together with Dain, prove [15] the existence of) a foliation of I + by special cuts (the 'nice 
cuts'), intended to model the 'system of rest frames'. Thus he is able to define not only a (supertranslation- 
ambiguity-frce) angular momentum, but higher order moments of the 'gravitational field', too. In their 
classic paper Bergmann and Thomson [21] raises the idea that while the energy- momentum of gravity 
is connected with the spacetime diffcomorphisms, the angular momentum should be connected with its 
internal 0(1, 3)-symmetry. Thus the angular momentum should be analogous with the spin. This idea was 
formulated mathematically by Bramson [22-24] , whose specific angular momentum expression was based on 
the superpotential derived from Hilbert's Lagrangian and the solutions of the asymptotic twistor equations. 
Recently, Katz and Lerer [25] could recover the Bondi-Sachs energy-momentum from the standard Noether 
analysis using asymptotically flat backgrounds, and suggested an expression for the angular momentum 
and center-of-mass. Interestingly enough, although these suggestions are based on different ideas and yield 
mathematically inequivalent definitions, all these resulting expressions bear some resemblance to each other. 
It is not clear whether the recent suggestion by Rizzi [26] , based on a special foliation of the spacetime near 
the null infinity and reducing to the ADM angular momentum at the spatial infinity [27], also takes such a 
form. 

At the quasi-local level the situation is even worse. Although the general framework of how the various 
quasi-local quantities should look like is more or less well understood (see e.g. [28]), only a few explicit 
definitions have been suggested for the energy-momentum and angular momentum. The first such suggestion, 
due to Penrose [12], is, however, not complete: In addition to the (symmetric) kinematical twistor A aj g two 
additional twistors, a Hermitian metric twistor h a| g> and the analog of the (skew, simple) infinity twistor lap, 
would be needed to reduce the ten complex components of h a p (the quasi-local quantities) to ten real ones, 
and to isolate four of them as the energy-momentum and six of them as the angular momentum. Although 
at the future null infinity these extra twistors exist and the construction works properly, reproducing the 
Bondi-Sachs energy-momentum and gives a definition for the angular momentum, it is not clear whether 
the additional twistors exist and the construction is viable quasi-locally or not. In principle the Brown- 
York approach also can yield both energy-momentum and angular momentum [29,30]. However, instead 
of a 4-covariant energy-momentum this approach yields separated energy and linear momentum, and it is 
not a priori clear whether they form a single Lorentz-covariant object. Furthermore, in lack of a unique 
prescription how the reference configuration should be find (imbedding of the 2-surface into a 3-planc [29] , 
or into the light cone in Minkowski spacetime [31,32], or restrict the imbedding in some other way [33]), this 
construction is not complete cither. To define Lorentz-covariant energy-momentum and angular momentum 
Ludvigsen and Vickers [34] used the Nester-Witten and the Bramson superpotentials, respectively, but their 
choice for the two spinor fields in the superpotential depends on the asymptotic structure of the spacetime. 
Thus their definition is not genuinely quasi-local. Genuinely quasi-local, manifestly Lorentz-covariant energy 
momenta, based on the Nester-Witten 2-form and either the holomorphic or the anti-holomorphic spinor 
fields, were suggested by Dougan and Mason [35], but they did not give any specific definition for the 
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angular momentum. Recently we suggested to complete the Dougan-Mason energy-momenta by a spin- 
angular momentum expression based on Bramson's superpotential, but, instead of the Ludvigsen-Vickers 
prescription for the spinor fields, we suggested to use holomorphic or anti-holomorphic spinor fields [36]. 
These have already been studied in various situations (pp- waves [36] and small spheres [37]). In the present 
paper these spin-angular momentum expressions will be calculated for large spheres near the future null 
infinity. 

In [22] Bramson introduced his superpotential as the superpotential for the conserved 0(1, 3)-current 
C° b , derived from Hilbert's second order Lagrangian (considering that as a function of the tetrad field 
and the 0(1, 3)-connection 1-forms). However, he defined C£ h as the partial derivative of Hilbert's action 
with respect to the connection 1-forms, while in gauge theories the conserved currents corresponding to the 
internal gauge invariance are the variational derivatives of the particle action with respect to the connection 
1-forms. In fact, this variational derivative is zero. Thus, in Section 2, first we show that Bramson's 
superpotential can be derived from Miller's tetrad action in a correct way, and then we review the status 
of our specific quasi-local spin-angular momentum expressions and provide the general formulae for the 
large sphere calculations. In Section 3. we review those structures of the future null infinity that we need, 
especially the BMS translations and rotations and their spinor form, and quote the asymptotic solution 
of the Einstein-Maxwell equations from [38]. Then, in Section 4, we calculate the anti-holomorphic spin- 
angular momentum for large spheres expanding that as a series of r" 1 . We will see that although this is 
diverging linearly and its finite order part is ambiguous in general, the Pauli-Lubanski spin vector built 
from the anti-holomorphic Dougan-Mason energy-momentum and spin-angular momentum is always finite 
and well defined. To demonstrate this, we need the 0(r _1 ) accurate expansion of the Dougan-Mason 
energy-momentum. Since, however, this calculation has been done only for stationary spacetimes [39] (and 
apparently an r _1 order term was overlooked), we have to clarify the asymptotic behaviour of the energy- 
momentum, too. In Section 5 we show that in stationary spacetimes the anti-holomorphic construction 
reduces to the 'standard' expression. The holomorphic construction, calculated in Section 6, is diverging 
quadratically in presence of outgoing gravitational radiation, but gives the same 'standard' expression for 
stationary spacetimes. Finally, in the Appendix we discuss how the various spinor equations determine 
special spin frames on two-spheres, in particular, on round spheres and on large spheres near the future null 
infinity. We show how the BMS translations and rotations can be recovered from them. We found that, 
in addition to the usual representations, they can also be recovered from the solutions of the limit of the 
Dirac-Witten equations on the unit sphere cuts of 1 + , too. We treat the 2-surface twistor equations also in 
the traditional way instead of the almost exclusively followed conformal method. 

Our conventions and notations are mostly those of [40]. In particular, we use the abstract index 
notations, and only the underlined and boldface indices take numerical values. The signature is -2, the 
Ricmann- and Ricci tensors and the curvature scalar are defined by — R a bcdX h Y c Z d := Vy{^zX a ) — 
Vz(VyX") — V[Y,z]X a , R a b '■= R e aeb and R := R a bg ab , respectively. Thus Einstein's equations take the 
form R ab - \Rg a b = -8nGT ab . 

2. Bramson's superpotential and quasi-local spin-angular momenta for spheres 

Let {E®,^*}, a = 0, ...,3, be an orthonormal dual frame field on an open domain U C M, and := 
V e -E£, the 0(1, 3)-connection 1-form on U. Then M0ller's Lagrangian [41] takes the form C := \Z\g\(R 
-2V a (i^ ab £ c c 7 c c b )) - j±c^\g~\(ElE» - i^hrf, where n ab := diag(l, -1, -1, -1). Considering the 
M0ller action as a functional both of the tetrad field and the connection 1-forms independently, Iu[E£, "Jabl := 
JjjL^x, for the 0(1,3) current we obtain 
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where a is some normalization constant, to be determined in some special situation. C° b is identically con- 
served, and we call the corresponding superpotential 2-form, := j^y^^c]' tne Bramson superpotential. 
The flux integral of this current for some spacclikc submanifold £ can be represented by the 2-surface integrals 
of the Bramson superpotential on its boundary S := <9E: J ab := -a § s W^\e ah cd = -ai£ ab c d f 5 W^f. 
If {E®} is a constant orthonormal (i.e. a Cartesian) frame field in Minkowski spacetime, then J ab is zero. 
To see this, it is enough to recall that in this case $ a = V a x a , the gradient of the Cartesian coordinate 
functions, and hence is an exact 2-form, whose integral on any closed orientable 2-surface is vanishing. 
Therefore, J ab is a measure of how much the frame field {E®} is distorted relative to a Cartesian basis in 
flat spacetime, i.e. a measure of the non-integrability of {E®}. (The current \/\g\C^ := SIu/SE^ gives the 
tensorial energy-momentum expression t b a for the gravitational 'field' found in [42,43] (and see also [44]): 

CI = Vlt\.) 

If S is any closed orientable spacelike 2-surface in M and A A , A = 0, 1, is a pair of smooth spinor fields 
on S such that they form a normalized spin-frame, i.e. £ AB \^)^ = e AB , where £ AB is the antisymmetric 
Levi-Civita symbol, then $ a :~ &a.b>^a^a> i s an orthonormal 1-form field on S, where the standard 

SX(2, C) Pauli matrices. Since the Bramson superpotential depends on $ a algebraically, the integral J ab 
can be expressed in terms of the normalized spin frame fields defined only on S. J ab is independent of the 
extension of {A A } off the 2-surface S. Translating the tensor name indices of J ab into spinor name indices, 
and defining its anti-self-dual part by J AA ' BB ' = : £ AB J A ' B ' +e A ' B ' J AB , we find that 

jAB =^£ x ?^ £a ' b '- (2 - 2) 

This is precisely the 2-surface integral of the spinorial Bramson 2-form [22]. Therefore, it has a natural 
Lagrangian interpretation in the sense that it is the spinor form of the superpotential of the conserved 
current derived from Miller's tetrad Lagrangian. In the present paper by quasi-local energy-momentum we 
mean an expression which is based on the 2-surface integral of the Nester-Witten 2-form: 

PAB := 8^g£^' Vbb ' XA ~ X b>Vaa>\%). (2.3) 

For any pair of spinor fields A A this defines a Lorentzian 4-vector, and it is natural to choose the spinor 
fields in (2.2) and (2.3) to be the same. (A more detailed discussion of these issues, and, in particular, the 
connection of these concepts and the current C a will be given elsewhere [45].) The quasi-local Pauli-Lubanski 
spin vector will be defined in the standard way by 

W := ^AA'BB'cc<DD'P BB '(e c ' D ' J CD + £ CD J C ' D ') = i(P A B< J* A' - ^a'b^a) ■ (2.4) 

(If the quasi-local mass to, defined by to 2 :— £ab£a'B'-P AA P bb , is not zero, then the dimensionally correct 
definition of the Pauli-Lubanski spin is ^Saa.'-) To complete the construction of these quantities, however, 
the spin frame field A A must be specified on, and only on S. 

In the present paper we will assume that the spin frame is holomorphic, fh e \7 e X A = 0, or anti- 
holomorphic, m e V e A J 4 = 0. Here m a and to" are the standard complex null vectors tangent to S and 
normalized by m a m a — —1, by means of which the metric area-element on S is — iTOr a TO^j . In fact, one 
can show that in the generic case there are two holomorphic/anti- holomorphic spinor fields and they can be 
normalized, provided S is homeomorphic to S 2 . Thus P AB and J AB that we will study here are the two 
energy-momenta of Dougan and Mason [35] , and the two spin-angular momenta that we introduced in [36] , 
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respectively. Then the formers can also be written as J> AB = —j^l^ ^(p'Aj^A 5 / + pXf^Xf,' )m[ a fh b ^ where 
7 = —1 for the holomorphic, and 7 = +1 for the anti-holomorphic spinor fields. In Minkowski spacetimc, 
for example, the independent holomorphic spinors are anti-holomorphic too, and they are just the constant 
spinor fields restricted to S, when obviously P AB = and, as we saw above, J AB is also vanishing [36,37]. 
If e\ := {o a ,l a }, A = 0,1, is a standard GHP spin frame adapted to the 2-surface, ejf := {— la,oa} is 

. . A A A 

its dual basis and the spinor components with respect to this frame are defined by A^" —: £~a^A' then the 
condition of holomorphy can be written as <9'A A + <r'A A = and <3'A A + /?A A = 0; and the condition of 
anti-holomorphy is equivalent to <9A A + ctA a = and $A A + p'A A = 0. Thus boldface capital indices are 
referring to a basis in the space of solutions, while the underlined capital indices to the GHP spin frame. For 
example, for a round sphere [39] of radius r the two linearly independent anti-holomorphic and holomorphic 
spinor fields are given by (A. 2. 2) and (A. 3. 2) of the Appendix, respectively. Substituting them into (2.2) we 
get zero, as could be expected in a spherically symmetric system. 

It is known that the anti-holomorphic Dougan-Mason energy-momentum is a future directed non- 
spacelike vector if S is the boundary of some compact spacelike hypersurface S, the matter fields satisfy the 
dominant energy condition on S, and S is convex in the weak sense that p 1 > [35]. (For the holomorphic 
construction the analogous convexity condition is p < 0.) Furthermore, if the dominant energy condition 
holds on the whole domain of dependence D(E) of £ then P AB is vanishing if and only if D(S) is flat, 
and this is also equivalent to the existence of two constant spinor fields on S with respect to the covariant 
derivative A e (see the Appendix); and P AB is null if and only if D(E) has a pp-wave metric and the matter 
is pure radiation, which is also equivalent to the existence of one constant spinor field on S [46,47]. J AB has 
already been calculated for (axi-symmetric) S bounding a pp-w&ve Cauchy development and it was shown 
that the Pauli-Lubanski spin Saw is proportional to the (null) Pab' [36] . 

We have already calculated (2.2) for the Ludvigsen-Vickers-, the holomorphic- and the anti-holomorphic 
spinors for small spheres S r of radius r with respect to an observer t a at a point o £ M [37]. Considering 
this to be a function of the radius, J^" B , it can be expanded as a power series of r. The leading term is 
— ^■Trr^TAA' BB't AA t B E e BF £^£p^ for the Ludvigsen-Vickers and the holomorphic spinors, whilst that is 
its negative for the anti-holomorphic spinors. In vacuum the leading term is — j^t & Taa' bb'cc DD't AA t BB 
t cc t D E s DF £ E A ~£p^ in all cases, where T a bcd is the Bel-Robinson tensor and {£&, £ A } is the dual Cartesian 
spin frame at o; i.e. E® = er AA ^a^A' an d $a = °aa' &a^A' f° rm an orthonormal dual frame at o. Thus 
the 'pure gravitational field' itself does not seem to contribute to the spin-angular momentum in r 4 order. 
On the other hand, for the leading term in the similar expansion of the quasi-local (anti-self-dual) angular 
momentum of the matter fields in Minkowski spacetime we get ^t^Taa'bb 1 

t AA' t B'E £ BF S (A s B) ^ 

coefficient rt A E e AF £^ s £p- ) is just an average of the approximate (and in Minkowski spacetime the exact) 
anti-self-dual boost-rotation Killing vector that vanishes at o, where the average is taken on the unit sphere 
S. Therefore, identifying this a.s.d. angular momentum of the matter fields with the r order holomorphic 
J AB , we obtain a = — 1, but identifying with the anti-holomorphic J AB we get a = +1. Thus there are two 
equally reasonable conventions: Requiring the r 4 order terms to coincide with that of the angular momentum 
of the matter fields, whenever the sign a is present in the basic definition (2.1) and in the two cases the 
r 6 order leading terms of J AB in vacuum are just minus the others; or to fix a once and for all for both 
constructions in the same way (e.g. by a = —1 as we did in [36,37]). We will see that in the latter case 
a = 1 would be a slightly more natural choice because the anti-holomorphic rather than the holomorphic 
constructions, both for the energy-momentum and spin-angular momentum, fit to the structure of the future 
null infinity. However, keeping in mind these observations, we retain our previous conventions, simply to 
facilitate the comparison with the previous works. 

In the rest of the paper the 2-surface will be assumed to be a large sphere of radius r near the future null 
infinity in the sense defined in section 3 below, and the corresponding spin-angular momentum will be denoted 
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by J AB . Then, according to the general philosophy of the large sphere calculations [38,39], we expand the 
spinor components A A , the GHP spin coefficients and the area element dS r as series of - . The expansion 
coefficients as functions of the remaining coordinates depend on the actual construction (holomorphic or anti- 
holomorphic). If therefore A A =: A A W + ±A A « + ^A A ( 2 ) + ... and dS r =: r 2 (l + ±s« + ± S W + ...)dS, 
where dS is the area element on the unit sphere S, then (2.2) takes the form* 

JAB _ _i_ y | r 2 ( A A(0) A B(0) + A A(0) A B(0)^ + ^ A A(0) A B(1) + ^(1)^(0) + ^(0)^(1) + ^(1)^(0)) + 
+ ( A A(0) A B(2) + A A(1) A B(1) + A A(2) A B(0) + A A(0) A B(2) + ^(1)^(1) + ^(2)^(0)) + j x 

x(l + i s W + l s ( 2) + -) d ^ 
\ r r z / 

(2.5) 

Therefore the r — > oo limit of J AB exists iff the r 2 and r order terms integrate to zero. Similarly, 

P r AB ' - -^j {r(p'WA^°)Af ) +p(i)A A ( )A?'(°)) + 

+ (p'< 1 >(A£< >Ag'< 1 > + A A «A B /(°)) +p>V\£m*;^ +p( 1 )(A A (°)A B '« + 

+ a a «a b /( o >) + p ( 2 )a a ^ b >))+ 

+ ^p'^{\^\B'(2) + A A(1) A B'(1) + A A(2) A B'( 0)) + p ,(2)( A A(0) JB'(l) + ^(1)^3' (0)) + 

+ p'( 3 )A A (°)A B /(°) + p« (A A (°)A B /( 2 ) + A A «A B '« + A A ( 2 )Af/(°)) + 

+ P w (a a(o) a b ' (i) + + p ( 3 )a a (°)a b ;(°)) + ...} x 

x(l + i s W + l s ( 2 ) + ...)d5, 

(2.6) 

which has a finite limit iff its r order term integrates to zero. 
3. The geometry of large spheres 

We assume that the spacetime is future weakly asymptotically simple [40], and let (it, r, C, C) be the standard 
Bondi coordinate system in a neighbourhood of the future null infinity 1 + based on a cut <S of 1 + (the 
'origin') (see also [48,49]). However, following [38], we adapt the null tetrad (and the spin frame) to the 
foliation of the outgoing null hypersurfaces by the 2-spheres <S( Ujr ) := {u = const, r = const}, which are called 
the large spheres. (In particular, the origin So above can be thought of as the r — > oo limit of the large spheres 
6>(o, r)-) Explicitly, let oaoa' := l a := V a w, which is geodesic, and impose the condition Z°V a os = 0. Let 
n a = l a l a be the other future directed null normal to <S( u , r ), normalized by n a l a — 1. Then m a = o A l A and 
fh a — l a 5 a are tangent to S( u _ r ), and the conditions above fix the basis. With this choice for the tetrad we 
have the following restrictions on the GHP spin coefficients k = e = p—p = t+/3'—(3 = p'—p' — t'—0'+(3 = 0. 
Therefore, ]d = (d/dr), df = (5 — (p — q)/3 — qr) f and d'f = (5+(p—q)j3—pf)f for any scalar / of type (p, q), 

* We use two different notations for the expansion coefficients: (i.e. when the index k is between 

parentheses) denotes the coefficient of r~ k in the expansion, which may turn out later to be zero. On 
the other hand, as is usual in the relevant literature, / k will denote the (A; + l)th nonvanishing expansion 
coefficient of the function / = However, for a function / both f k and will never appear. 
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where 5 := m a V a =: P(d/d() + Q(d/dQ, and £, £ are the standard complex stereographic coordinates on the 
2-sphere. In particular, on the unit sphere in the Minkowski spacetime the edth and edth-prime operators 
take the form df := ^(1 + (()(df/dO + ^(p - q)Cf and d' f := ^(1 + QQ{df/dQ - ^{p - q)(f, 
which will be used in the subsequent calculations. Overdot will denote partial derivative with respect to u. 

In the coordinate system (u, w,(, 0, where w :— r" 1 and the future null infinity is given by w = 0, 
the BMS vector fields have the form k a = (H + (c 1 + c% u)(^) Q + c i 6(^)° + c i ^ i {-^) a + C^r" 1 ), where 
H is an arbitrary real valued function of £ and £, c 1 € C, i = 1, 2, 3, and t\ := — (( + £)(1 + (C)^ 1 , := 
i(C - C)(l + CC)-\ t 3 := (1 - CC)(1 + CCr 1 and & := (1 - ( 2 ), 6 := i(l + C 2 ) and & := K For the generators 
of the BMS supertranslations c 1 = 0, and, in particular, the independent BMS translations take the form 
*a = *a(^) a +0(r _1 ), a = 0, 3, where to := 1 and U for i = 1, 2, 3 are given explicitly above. The functions 
i a can be written as i a = er AB tatb>, where r := exp(ia) v^2(l + CC) ^ an d T i '■— ~ exp(ia)v / 2C(l + CC) - ^ 
and exp(ia) is an unspecified phase. In Minkowski spacetime the standard constant orthonormal frame field 
{E£} (i.e. the translational Killing vectors K£) has precisely this asymptotic form, and if the phase exp(ia) 
is chosen to be — i then the functions ta are just £^oa, the contractions of the Cartesian spin frame with 
the GHP spin vector oa- In fact, in terms of the GHP spin frame e A = {o A ,i A } the Cartesian spin frame 
£ A = {O a , I A } has the form 

Thus the BMS translations above are properly normalized. In the conformal approach (see, e.g. [40]) 
the GHP spin frame is rescaled by d A = Q~ 1 o A and i A = i A , thus expressing o A and i A by the non- 
physical spin vectors 6 and t in (3.1), for the asymptotic behaviour of the Cartesian spin frame (i.e. the 
spinor constituents of the translation Killing vectors) we directly obtain £ A = tj±l a + 0(0), where the 
exp(ia) = — i choice was made and the conformal factor is ft — r _1 . To find the proper normalization 
of the boost-rotation BMS generators too, let us consider the Minkowski spacetime again. The anti-self- 
dual part in the name indices of the boost-rotation Killing vectors K£ b of the Minkowski spacetime is 
K ab = 71 o'AB^oi + | £ i jk - K j°k)> whcre iJ> k = 1, 2, 3 and er A B := V2ct a B '(t b 'b, the standard SJ7(2)-Pauli 
matrices, and eyk is the alternating Levi-Civita symbol, whose indices are moved by the constant (negative 
definite) metric := diag (— 1, — 1, — 1). They yield on J+ the BMS vector fields fc AB := ^ AB li+ = 
^a AB (?xt i (| r ) a +e i (^) a ) + 0(r- 1 ) = ^(^ut^'+^l + CO^i^^ + Oin). Note that the functions 
6(1 + CC) 1 can be written as &(1 + £C)~ 1 — — ct AB TaTb, where the phase in r A has been (and in the rest of 
the paper will be) chosen as exp(ia) = — i. Therefore, apart from the (in general super) translation content, 
the generator of the anti-self-dual rotations at I + is 0AB£i(l + CC) _1 = -T( A t B ), i.e. minus the symmetrized 
product of the functions ta- In the Appendix we discuss how the generators of the translations and of the 
anti-self-dual rotations can be recovered from the solutions of various spinor equations on the cuts of X + . 

Suppose that, at least in a neighbourhood of X + , the only matter field that we have is the electromagnetic 
field, represented by the components of the Maxwell spinor (see e.g. [48,49]). Imposing slightly stronger fall- 
off conditions than that coming from the definition of the future weak asymptotic simplicity, viz. assuming 
that ^0 = r^ 5 ip a ° + r^ipo 1 + 0(r~ 7 ) and (j> = r~ 3 o ° + r -4 ^ 1 + 0(r~ 5 ), Shaw found the asymptotic 
solution of the Einstein-Maxwell equations [38]. We need this solution with accuracy 0(r~ 3 ). In particular, 



(3.2) 
(3.3) 
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This implies that the area element of the large sphere of radius r is d<S( U;r ) = r 2 (l — r 2 <t <t + 0(r 4 ))d<S; 
i.e. in equation (2.4) — and = —a°a . The spin coefficients with definite (p, q) type are 

a=4 7( r + O(r- 4 ), (3.4) 



P = -- r -^°°° + 0(r-% (3.5) 

+ cx ( dV + o^'V) + &a° &a° ~ 6a° 6a° G<h°fo°) + 0(r~ 4 ), (3.6) 
P'=h + h (^2° + + o^V) - 1 ( da° J a - \a°a° + a d' da°+ 

+ a d d'a° + \ o^Vi° + \ o^i'° - 2G^°^°) + 0(r~ 4 ), (3.7) 
r =^od'<7° - ^(2a\da° + +0(r- 4 ); (3.8) 
while the spin coefficient representing the GHP connection 1-form is 

Finally, for the Weyl and Maxwell spinor components one has the familiar peeling: ip n = r n ~ 5 tp n ° + 0(r™~ 6 ), 
n = 0,...,4, and = r"- 3 0„° + C>(r"- 4 ), n = 0,1,2, where V>4° = "03° = - and V 2 ° - V^ = 

<7°a - a ,?' + d' 2 a° - d 2 a°. 

The spacetime will be called stationary, if it admits a timelike Killing vector field, at least on an open 
neighbourhood of I + , which can be extended to X + into a BMS translation [48,38]. Then there is a cut, 
chosen to be the new origin, whose asymptotic shear is vanishing, and the asymptotic shear on the cut given 
by the supertranslation S : S — > R with respect to the new origin is a = — d 2 S. Then, after an appropriate 
translation of the new origin, the asymptotic value of the Weyl and Maxwell spinor components that we 
need in what follows on such a cut take the form ipi° = — 3G od(MS + U), tp2° = —GM and <f>i° = ^(e + i/z). 
Here M, e and \i are real constants, interpreted as the total mass, the total electric charge and the total 
magnetic charge, respectively, and J is a real function with structure J = X)m=-i JmXi.m for some constants 
Jo and J±i, where Y\, m are the standard j = 1 spherical harmonics. Rewriting J by the familiar polar 
coordinates (9,(f>) (defined by £ =: cot | exp(i^)) into the form J = ji sin 6> cos </> + 32 sin 9 sin <j> + j'3 cos 9, 
the real constants j\ — {31,32,33) defined in this way are interpreted as the components of the total spatial 
angular momentum of the stationary solution. In particular, for the Kerr-Newman solution j\ = 32 = 
and j 3 = Ma, and hence on the shear-free u = const, cuts = — ^iGMasin#exp(i<^). (Its apparent 
deviation from the expression given in [17,38] by the factor — exp(i</>) is a consequence of the different choices 
for the holomorphic coordinates on S: Our choice is ( [even if spherical polar coordinates are present], while 
in [17,38] it is -logC = log tan f +i<f>.) 



4. The anti-holomorphic spin-angular momentum 

Substituting the expansion =: Xa (°) + ^A^ W + -^Xa ^ +••• of the spinor components and the expressions 
(3.2-9) for the functions P and Q and the spin coefficients into the equations defining the anti-holomorphic 
spinor fields on S( u , r ), we obtain the following hierarchy of equations 
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o 0Ai(°> + ^Ao< o >=O, (4.1.a) 
o#A (0) = 0, (4.1.6) 
dX^ + h ( V = -(V2 + t^' + o^a°)A ( \ (4.2.o) 
o$A (1) - 0, (4.2.6) 

o^ 2 ) + iAo^ = a d'X^ - (^ 2 ° + <7°a'° + dV)A « + (a + A 1 <°) + 

+ (o&7° o^V° + a° + a o^o^'a + \ + \ Mv* - 2G<h°fo°) X {0 \ (4.3.a) 

dA ( 2 > = ^(o^'Ao^) - A^ 1 )) - (a\da° + ^! )A (°). (4.3.6) 

Thus by (4.1.a-b) the zeroth order spinor components are just the components of the constant spinors of 
Minkowski spacetime, and hence, in addition to (4.1.a-b), they satisfy od'Xi^ = and q9'\o^ = Ai^ ^, too. 
Therefore, there are precisely two solutions A A W), A = 0, 1, of (4.1.a-b), and we choose them to be given 
explicitly by (A. 2. 2) with p' = ^- (see the Appendix). 

Since the left hand side of (4.1), (4. .2) and (4.3) are the same, the solution in each order will be the sum 
of a particular solution and the general zeroth order (i.e. the constant) solution. Therefore, the general r~ 2 
accurate anti-holomorphic spinor fields form a six rather than the expected two complex dimensional space. 
In fact, in the small sphere calculations [37] we had similar spurious solutions, due to the fact that there is no 
natural isomorphism between the space of the anti-holomorphic spinor fields on two different two-surfaces. 
Thus they represent a gauge ambiguity in the first and second order corrections, and the physical quantities 
must be invariant with respect to the substitutions A^ W A^ ^ +CaA a ^ and A^ ( 2 ) ^ Xa ^ +D a \^ ^ 
for any constants Ca and Da- 

To compute the spin-angular momentum based on equation (2.5), first observe that the integral of the r 2 
order term is vanishing, because that is just the spin-angular momentum in Minkowski spacetime. However, 
the r order term of the integrand is 

A^AfW + A^WAfW + A A Wa^ 1 ) + \?W)$<n = 

= -2 3(Xt (0) Af « + \fM\?W) + 2A A (°) ( od\fW + iA^D) + 2A? (°) ( d\f « + ^A^ 1 )) = 
= -2od(\? (°>Af « + Xf^XfV) - 2(V + <j°<j" + d 2 * ) (A A (°)A?(°) + A?< >A£< >), 

where we used equations (4.1) and (4.2); furthermore, by (4.1) the last of the second term is a total divergence: 



6 2 a°(\?W\*M +A B(0 )A A(0)) = 5( ^0 (A A(0) A B(0) + ^(0)^(0)) + -0 A A(0) A B(0)y (<u) 

Thus for the anti-holomorphic spinor fields J AB is diverging at 1 + unless the integral 

L AB :^--L^(^ 2 + (T 0a0)(Af( )A B ( )+A A ( )Af( ))d5 (4.5) 
is vanishing. (We return to the discussion of L AB below.) Next let us consider the r° order term of (2.5): 
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\ j (A^^Af < 2 ) + A^Af^ + A a < 2 )A? (°) + A A (°>A^ 2 ) + Xf W^W + Af^ )- 

- ^(A^Af W + A A (°)A?(°)))d5 = 

= £( - (o^Af (°))A? < 2 > + Af«(- c0A A « (^ ° + ^ + 3V)A£(°)) + ^A A ( 2 )Af (°)- 

- ( o^Af A A < 2 > + A A « (- odAf« - (^o° + <r°d° + d 2 a°)X^) + Ia?< 2 >A^°>- 

_ 1 CT 0-0( A A(0) A B (0) + A A(0) A B(0)^ d<S = 

= £(Af(°)( 5Af ( 2 ) + 1a b < 2 ) - ±aVA?<°>) + A? (°VAf < 2 > + Ia a < 2 > - i^A^)- 

- (^o° + a ^ + d 2 a°)(\^\f + A A «A B (°)))dS = 
= £ ( - //< 3 >(A A (°>A? (°) + Af (°)A B (°)) + 

+ (Vi° + 2,7° &7°)A A (°)AfW + <7°(A A (°> d'A B « + Af (°) tf\? «)- 

- (^2° + <7°<7° + dV)(A^°)A? C 1 ) + AA(1) A B(0) + A A(0) A B(1) + ^(1)^(0))^ (46) 

where p'^ is the 3rd order term of p' in (3.7), and we used (4.1. a), (4. 2. a) and (4. 3. a). Substituting (3.7) 
here, using (4.1.b) and the consequences d'A A (0) = and 5'A A(0) = A A W of (4.1), the integral of the first 
three terms of the right hand side of (4.6) can be written as 

j ((Vi° + 2 ( 7 t Va°)A A (°>A? (0) - ^V°(A A (°)Af« + Af ^Af «) + ( da\d'a - \<J°a + 

+ a d' da° + a d 8'a° + \ o^V + \ o^/ - 2G<h.°ti>°) (X^Xf^ + A A (°)A B (°)))dS = 

= j ((1^,0 + ^,,^0^(0)^(0) _ o5 ' CT 0( A A(0) A B(l) + A B(0) A A(1))_ 

- ( #*° da° + ^"a + 2G0!%'°) (A A (°>A B (°) + A A (°)A B (°)))d5 = 

= \ jf ((^° + 2<j" od'*° + d'(a°a ))X^X^ AG J '(^ )(A*<°> A b W + A A (°>A B (°))- 

- 2 £>V°(A A (°)(A B « + A B (°) ^ ) + A B (°)(A A (D + A A <°) 0*°)))dS. (4.7) 

Using (4. 2. a) again, the integral of the last term on the right hand side of (4.6) can be rewritten as 

- <f (^ 2 ° + o°d*> + od 2 v°) (a a (°)A b « + A A «A B (°) + A A (°)A B « + A A «A B (°>)dS = 
•J s 

= £(A B «( dA A « + ^A A «) + A A ( 1 )( 5Af « + iA B «)- 

- (^ 2 ° + a°a0 + d 2 a°) (A A «A?(°) + X^X^))dS = 
= £{A A «(iA B ( 1 )-(^ + a°aO+ /s°)Af(»)) + 

+ A B(1) (V>2° + <r°*° + odV)A A (°)) }dS. (4.8) 
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Therefore, A A ' 2 -* and A A ( 2 -* are not needed to calculate the anti-holomorphic spin-angular momentum, but 
X A ^ and A A ' 1 - ) do appear in (4.6) explicitly. Since, however, the physical quantities should not be sensitive 
to the addition of the spurious zeroth order solutions to A A ^ and X A ^\ and by such solutions A A ^ = can 
always be achieved, whenever (4.8) gives zero, that should be a gauge term. To see that, in fact, this is the 
case, recall that (4.1.b) and (4.2.b) are the same, thus we may write A A W = C A C A^°) = -2C A co dXf (°) 
for some constant complex 2x2 matrix C A c- (However, C a b is not quite arbitrary, that is restricted by 
the requirement that the pair of anti-holomorphic spinor fields should form a normalized spin frame. In 
fact, from e AB = (A A (°) + ±A A « + ...)(A|(°> + ±A|« + ...) it follows that C AB must be symmetric.) 
Substituting this into (4.8) and using (4. 2. a) we obtain 

-2(C A C 5E + C B C 5£) j odAfW^A?' 1 ) - (^ 2 ° + a°a° + *V) A? '(°>) dS = 
= (C A C S B + C* C 5 A ) j Ap(°)(o5A?« - 2 3(V2° +^'° + o0V)A?<°> + 

= (C A c <5g + C B C <5 A ) £a?W (-^A?^) - 2 % 2 ° + <7°*° + o 2 «7°)A?<°))dS = 
= \(C A C C B » + C B C C A D ) £ A C (°)A? (°)d5+ 
-(C A C 5 B + C B c 5 A )£(ip 2 ° + o°d° + d 2 a°) (a c (°)A 1 d (°) + Ai c < >A D < >)dS. (4.9) 

However, the first integral on the right hand side of (4.9) is vanishing, while, taking into account (4.4), the 
second is seen to be proportional to L AB above. Therefore, the anti-holomorphic spin-angular momentum 
for the large sphere of radius r is 

J AB = rL AB + (C A C 5* + C B C 5 A )L CD - E AB + / AB + 0{r~ l ), (4.10) 
where we introduced the notations 




(4.11) 



/AB := 8^G jf {(&'° + 2*°o0V+ d'(a°a ))x A ^X^- 

- 2 d'a° (A A W (Af « + A»<°> da°) + A?<°> (A A(1) + A A W da°))}dS. (4.12) 

Therefore, as we noted above, the quasi-local spin-angular momentum based on Bramson's superpotential and 
the anti-holomorphic spinors is diverging at the future null infinity, furthermore, its finite part is ambiguous 
unless L AB is vanishing. In addition, contrary to expectations, the electromagnetic field contributes to 
J AB in the 0(1) order. On the other hand, the first three terms together in the integrand of I AB is just 
the integrand of Bramson's specific spin-angular momentum expression based on the asymptotic twistor 
equation [22] . Though in the present case the spinor fields A A are anti-holomorphic and not the solutions of 
the asymptotic twistor equations, the zeroth order part of A A , which appears as the coefficient of the first 
three terms of the integrand of / AB ; coincides with the zeroth order part of the solutions of the asymptotic 
twistor equations. Furthermore, the coefficient A A (°)A B (°) is just minus the generator of the anti-self-dual 
rotations on the u = cut (the 'origin') of 1 + (see Section 3. and the Appendix). It might be worth 
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emphasizing that these generators of the anti-self-dual rotation BMS vector fields emerged naturally, like 
the approximate rotation-boost Killing vectors in the small sphere calculations [37], without putting them 
into the general formulae by hand. Although J AB depends on the solution A^ 1 ) of (4. 2. a) (and hence its 
integrand is a genuinely non-local expression), that is independent of the gauge solutions. In fact, using 
o5'A A (°) = 0, it is easy to see that the addition of a gauge solution to A A ^ changes the integrand by a total 
o^'-derivative. To clarify the meaning of the term of the integrand involving A A W, recall that the BMS 
vector fields k^u are tangent only to the origin cut, and they can be represented completely by A A (°)A B (°) 
only there. On general cuts fc^ B contain correction terms proportional to the generator {-§^) a of 2 + . Then 
the integral of the term involving A^*- 1 ) can be considered as the contribution of this correction term built 
in a non-local way from the fields on the actual cut and the spinor constituents of the BMS rotations on the 
origin cut. 

However, L AB is not zero in general, because its components are proportional to that of the spatial part 
of the Bondi-Sachs energy-momentum 

ooP AB ' :=- i ^£(^ + a aO)A A ( )A B '(°)d5. (4.13) 

To see this, let us substitute the explicit solutions (A. 2. 2) into (4.5). We obtain L 00 = -\/2oo-P 01 \ L 01 = 
-^(oo-P 00 — oo-P 11 ) an d L 11 = V^oc-P 10 , i.e. L AB represents the linear momentum. In fact, for A = 
B = 0, A = B = 1 and A = 0, B = 1 the coefficient of t/> 2 + in (4.5) is -2£(1 + CC)~\ 2<C(1 + CO -1 
and (CC — + CO -1 ; respectively, which are proportional to the independent spatial BMS translations. 
Therefore, the anti-holomorphic spin-angular momentum can be finite only in the center-of-mass system 
(i.e. when the spatial components of the Bondi-Sachs energy-momentum are vanishing), and hence / AB in 
the 0(1) part of (4.10) appears to represent only the intrinsic or spin part of the total angular momentum, 
while ri AB appears to be the orbital part of the angular momentum. To check whether this interpretation 
is correct we should calculate the quasi-local Pauli-Lubanski spin (2.4) built from the quasi-local anti- 
holomorphic Dougan-Mason energy-momentum P AB and the anti-holomorphic spin-angular momentum. 
However, to compute the spin, we need to know the Dougan-Mason energy-momentum for large spheres with 
0(r _1 ) accuracy. Since this has been calculated only in stationary spacetimes [39] (where a physical term 
was apparently overlooked and the gauge ambiguity caused by the spurious solutions was not considered at 
all), first we must calculate this. 

The r order part of (2.6) is vanishing, because $ s A A (°) Ag' (°)dS = 2 $ s A A (°)A B '(°)dS = 47rcr AB ' (sec 
Appendix A. 2), i.e. P AB has a finite r — > oo limit at T + . Substituting (3.5) and (3.7) into the finite term of 
(2.6), using (4.1.a), $'A A(0) = A A <°) and its complex conjugate, (4.2. a) and the fact that ip 2 ° + <7°&° + d 2 a° 
is real, we obtain 
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Since the last term of the integrand is a total d-dcrivative: ( dV)A A (°) A B ,'(°) = d( a da° A A (°) A B / (°)) - 
(o5a°)A A (°) dA B /(°) = d( da° A A (°)A B '(°) - a°X A ^Xf,'^), the finite part of (2.6) is, in fact, the Bondi- 
Sachs energy-momentum (4.13). 
The Oir- 1 ) term of (2.6) is 

£((p /(3) - ^°a )A A (°)A B '( ) + + <j°a° + d 2 a )(X A ^\B'(i) + A A(i) X B'(o )) + 

+ ^(A A(0) A B /( 2 ) + A^A^ 1 ) + A A ( 2 )A B '(°)) - 5'A A(0) A B '( 2 ) - X A ^ d\%' w - A A « A B '«)dS = 
= jf ((p'^ - ±a°a )X A m*'W + + a o- + o ^_ 0) (a a ( o) X b/(i) + a a ( i) X b/(o )) + 

+A A (°)(o5'A B '( 2 ) + 1a b '( 2 )) + A B '(°VA A ( 2 ) + ±X A ^) + iA A «A B '« - A A ( 1 )A B '( 1 ))d5 = 
= j> (-{2G<h°fa> + oBa° od'a°)X A m*'^- 

-{ da°)X A m?;W _ ( ^V°)A A ( 1 )A B '(°) + ^A A(1) A B ' (1) - A A «A B '«)dS = 
= I (-2G^ ^°X A ^XV'(o) + iA A «A B '« - (A A « + A A W 3a°) (Xf,' {1) + A B '(°> d'a°))dS. 

(4.14) 

Here first we used $'A A(0) = A A (°\ and then (4.3. a), (4.1.b), (3.7) and the fact that V> 2 ° + a°& a + o^V 
is real. However, as we saw above, X A ^ — C A cA^°) = — 2C A c o^Af (°\ and by means of which (4. 2. a) 
takes the form d(X A ^ + X A ^ da°) = C A C dXf^ - (^ 2 ° + a°a )X A ^\ Since d acting on s = -\ 
spin weight quantities is isomorphism (see e.g. [47,50]), we can write A A ^ + A A (°) da-° = C A cAp^ ^ — 
O0 _1 (W>2° + a°a a )X A ^). Substituting these into (4.12), using d'X A ^ = A A (°), (4.13) and the expression 
for the L 2 -scalar product of the components A A (°) and A A (°' given in Appendix A. 2, finally we obtain 

P r AB ' = ocP AB ' + i(C A c<5g: +C B d^^oo P CD ' - ±F AB ' + Im ab ' +0(1), (4.15) 
where we used the notations 

^ AB ' £^°^' A A(0) ^' (0) d5, (4.16) 
M AB ' := - ^ £ otT 1 ((^° + a°>)A A (°)) o^'" 1 ^ + A B '(°>)dS. (4.16) 

Therefore, as it could be expected, the r _1 order term in the expansion of the anti-holomorphic Dougan- 
Mason energy-momentum is ambiguous, and, in addition to the electromagnetic contribution, there is an 
extra term. Note that E AB above is related to F AB ' like L AB to oo-P AB ', and hence E 00 = -\/2F 01 ', 
£0i = -i= (F 00 ' - F 11 ') and E 11 = V2F 10 ' . 

Substituting (4.10) and (4.15) into (2.4) and using how the components of L AB and E AB are related 
to those of oo-P AB and F AB , respectively, we obtain 

Sab' = i( M ^/ C B- + M AC 'L C b' - ocPb'c/ C a - M B , c L c A y (4.17) 

The diverging term, the ambiguities and the contribution of the electromagnetic field disappeared. Therefore, 
the quasi-local Pauli-Lubanski spin vector built from the anti-holomorphic Dougan-Mason energy-momentum 
and the anti-holomorphic spin-angular momentum (2.2) has finite limit at the future null infinity, it is free 
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of ambiguities, and is built only from the gravitational data. Note that our Pauli-Lubanski spin is free 
of the so-called supcrtranslation ambiguity, because this is defined in terms of the solutions of an elliptic 
differential equation on the cut in question, and not by means of the BMS boost-rotation vector fields. Thus 
the present construction is similar in its spirit to that of Penrose [12]. The fact that we could derive only 
a Pauli-Lubanski spin is compatible with the idea of Bergmann and Thomson [21] that the gravitational 
angular momentum should be analogous to spin (justifying the 'spin-angular momentum' terminology), but 
raises the question as whether that should be completed by an orbital angular momentum part or not. 
Another interesting issue is how the Pauli-Lubanski spin changes under (infinitesimal) supertranslations, in 
particular, under time translations; or whether there exists a flux for Saa' through 1 + or not. However, 
these questions are beyond the scope of the present paper. 



5. Stationary spacetimes 

Suppose that the spacetime is stationary. First we show that all the terms of the integrand of (4.12) involving 
the asymptotic shear together integrates to zero. Bramson already showed that 2a° d'a° + od'(a°a°) — 
d' 3 {S d 2 S + ±( dS) 2 ) - d(S a d' 3 dS + 3 d' 2 dS dS + 3( d'S) 2 ), which, together with dX A ^ = 0, 
od' A A (°) = A^ -* and q3 A A (°) = 0, gives that the first two <7°-terms of the integrand of (4.12) give zero. To 
evaluate the last term of the integrand, first we must solve (4. 2. a). In stationary spacetimes (4. 2. a) takes 
the form dA A « + iC A B A^ ) = (GM + Q d 2 Q d' 2 S)X A< -°l Using dimker<9 ( _ l!0) = 0, this can be solved 
explicitly Its solution is A A W = „0( Q d' 2 S)X A< -^ + (-2GMS A + C A B )A B (°). Substituting this into the last 
term of the integrand, that takes the form 

- 2 J a (a a (°) (Af « + A*W d<T°) + Af (°> (A A W + X A ^ = 

= - 2((-2GM5£ + C A C )S* + {-2GMbt + C B D )5 A ) o 0Va?( o >A?<°>). 

Therefore, the last term of the integrand of (4.12) integrates to zero, too. By -02° = —Gm the integrand 
of i AB is the total (^'-derivative d'(-GmX A ^ A B W), and by fa = \{e + ut) the integrand of E AB is 
d'(j[e 2 + ^ 2 ]A A (°)Ajf (°)). Thus in stationary spacetimes (4.10) reduces to 

J r AB = £^ A A (°)A B (°)d < S + O(r- 1 ), (5.1) 

which is the 'standard' expression for the angular momentum in stationary spacetimes. In fact, substituting 
V>i° = -3G d(MS + i£m=ii J m Y hm ) here and using & A A (°) = A A (°\ wc find 



tAB _ 

r ~ 8ir 



3M £s(x A ^Xf^+X A ^X^)dS-^ £ X:£^(A A (°)AF( )+A A ( )A B (°))d5 + O(r- 1 ). 



m— — 1 



Substituting the specific solutions (A. 2. 2) here and defining c m := M J> 5 Yi. m SdS and the corresponding real 
components c 1 = (c 1 , c 2 , c 3 ) by Y^,m=-i c ™^i,m —■ c 1 sin 8 cos (f> + c 2 sin 9 sin 4> + c 3 cos 9, for the leading terms 
we finally obtain J lj = e^ju and J ok = — c k . Thus J AB reproduces the jk's, and the c k 's can be interpreted 
as the components of the relativistic center-of-mass. In particular, for the Kerr-Newman solution on the 
shear-free cuts J 12 = —Ma. Using the expressions of Appendix A. 2 for the L2 scalar product of the spinor 
components A A (°\ we obtain ooP AB ' = Mcr AB ', F AB ' = ±(e 2 + M 2 ) cr AB ' and M AB ' = -2GM 2 a AB ' . 
Hence in stationary spacetimes the Pauli-Lubanski spin reduces to that of Bramson. As he showed [23], this 
is invariant with respect to supertranslations. 
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6. The holomorphic spin-angular momentum 

For the components Xa =■ Xa ^ + ^Xa ^ + 77 Aa ^ + ••• of the holomorphic spinor fields on <S(„. r ) we obtain 
the equations 



oc?V 0) -Xi (0) -0, (6.1.a) 

o^'Ai^ -a A (°) =0, (6.1.6) 

od'Xo^-X^ = o%°A (0) ), (6.2.a) 

od'X.W _^o Ao (D =a°( BX 1 W + ^A (°)) - ^((o^°)Ao (0) ), (6.2.6) 

tfV 2) - A/ 2 ) = %°Ao W ) - + <x o^^°) A (0) , (6.3.a) 

ocJ'A^ 2 ) - a"A ( 2 ) = ^(ocJA^ 1 ) + 1a «) - ^((o^°) A (1) ) + + a dV) X^ + 

+*° od*° dX ^ + (IaV 2 + ^o^i'°+ o^(^o^°) +a%c5V-G0 2 o W o )A o (°). (6.3.6) 

Therefore, the zeroth order holomorphic spinors are noi constant unless &° is vanishing, and hence the coef- 
ficient of the r 2 order term of (2.5) is not zero in general. Therefore, the quasi-local spin-angular momentum 
(2.2) based the holomorphic spinor fields is diverging near2 + in presence of outgoing gravitational radiation. 
Thus let us concentrate on spacetimes for which a = 0. Then Xq^°\ Ai^ ^ are components of a constant 
spinor field on S, implying the vanishing of the r 2 order term in (2.5), and the vanishing of the first term on 
the right hand side of (6.2.b). Then, however, its general solution is A A W = -( <fc°)A A(0) + C A B A B (0) for 
some 2x2 complex matrix C b, by means of which the solution of (6. 2. a) is A A W = C A B A B (°). Therefore, 
the integrand of the r order term of (2.5) is 

A A <°)Af W + A A «Af (°) + A A (°>A B « + A A «A B (°) = 
= (5 A C B D + C A c^)(A c ( >A?(°) +A D (°)Af(°)) -2o^( ( t A a (°)A b (°)). 

However, its integral is zero, i.e. in absence of outgoing gravitational radiation (i.e. if & — 0) the quasi-local 
angular momentum based on (2.2) and the holomorphic spinor fields has a finite limit at the future null 
infinity. 

To calculate this finite value let us consider the r° order term of (2.5). Using the explicit solutions for 
A A W and A A< ^ above, the fact that A A (°) is constant and (6. 3. a) we obtain 



£( 



A A(0) A B( 2) + A A(1) A B(1) + A A(2) A B(0) + A A(0) A B(2) + ^(1)^(1) + A A(2) A B(0)_ 

- a°a (A A ( )A B (°) + A A (°)A B (°)))dS = 

/((o5'A B ( ))A A ( 2 ) + (o5'A A (°))A A ( 2 ) + 

+ C A C A^°) (- ^°A B (°) + C B D A?(°)) + C B D A?(°) (- ^°A A (°) + C A c Af < >) + 
+ A A <°)A B ( 2 ) + A B (°)A A ( 2 ) - a°a°(A A (°) Q d'X^ + A B (°> Q d' X^))dS = 

A A(0) (- o^'A B(2) + A B < 2 )) + A B W (- o$'A A(0) + Af( 2 )) - a°a° £>(A A(0) A B < )) - 

- (C A c <5g + C B D <5 A ) o^°A c ( )A D (°) + C A cC B D (A c (°)A?(°) + X^X^))dS 
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which is precisely Bramson's specific spin-angular momentum expression based on the asymptotic twistor 
equation. Thus, substituting cr° = — od 2 S here, finally we obtain (5.1), i.e. for stationary spacetimes the 
holomorphic and the anti-holomorphic constructions give the same 'standard' expression. We expect that 
at the past null infinity the holomorphic construction works properly, and the anti-holomorphic diverges in 
presence of incoming gravitational radiation. 



Appendix: Special spin frames and asymptotic symmetries of J+ 

On the bundle S A (<S) of unprimed spinors over S two connections can be introduced in a natural way by means 
of the spacctimc connection. The first is the simple projection to S of the spacetime covariant derivative 
operator, A e := E^V/, where the orthogonal projection to S can be given e.g. by iljj = —m a fhi, — fh a mb- 
The other is the spinor form of the covariant derivative operator 5 e defined on any spacetime vector field X a 
by S e X a := Yl^A e (U b c X c ) + (6% - ng)A e ((5* - n b c )X c ). The difference of these two connections is essentially 
the extrinsic curvature of S in the spacetime. (For the details see e.g. [50,51].) 

Representing the spinor bundle by the line bundles E(p,q), p — q e Z, the derivative operator S e 
can be represented by the edth and edth-prime operators d and d' [52]. Explicitly, by the isomorphism 
Sa(<S) w E(l, 0) © E(— 1, 0), A^ w (Ao, Ai), the derivatives m e 8 e \A and fh e S e XA are represented by the cross 
sections (5Ao,$Ai) and {&' Xq, d'Xi), respectively. Denoting e.g. 3 acting on the cross sections of E(p,q) 
by d( P ,q), on topological 2-spheres dim kerf ( P!P+ „) = dim ker d'( p + n ,p) = f° r an y P ^ R and n e N, and 
dimker<9( p+ „.p) = dim ker d, pp+n \ = 1 + 2n for any pgR and n = 0, 1,2, .... (For the corresponding kernels 
on higher genus surfaces see [50].) The irreducible parts of A e and S e , respectively, are 



T-(Ao.Ai) 
A+(A ,A!) 
-A-(Ao.Ai) 

-r+(\ ,Xi) 



--m e (A e X A ) o A = 3X + aXi , 
=?fi e (A e A A )o /l = d'Xo + pX u 
--m e (^ e X A )i A =dX 1 + P 'X^ 
--m e (A e X A )i A =d'X 1 +a'X Q , 



i"(A ,Ai) := m e (6 e X A )o A = $A , (A0.1.a,6) 

S+(X , Ai) := rh e (5 e X A )o A = &'X , {A.0.2.a, b) 

-<T(A ,Ai) := m e (5 e X A )t A = (A0.3.a,6) 

-* + (Ao,Ai) :=fh e (6 e X A )t A =3'X 1 . (A0.4.a,6) 



In particular, a spinor field A^ is constant on S with respect to A e iff (Ao, Ai) G kcr(T _ © A+ © A~ © T+). 
As Bramson showed [53] , such a spinor field does not exist even on large spheres near the future or past null 
infinity in a general asymptotically flat spacetime; and for a finite topological 2-sphere, being the boundary 
of some compact spacelike hypersurface S on which the dominant energy condition is satisfied, the existence 
of a constant spinor field implies that £)(£) must have a pp-w&ve metric [46,47]. To weaken the notion 
of the A e -constant spinor fields on S or on the cuts of 1 + , in principle there are six natural possibilities. 
These are represented by the kernel of the first order operators (see [47]) A := A+ © A~, H~ := A~ © T~, 
7i + := A+ © T+, C+ := A+ © T~, C~ := A~ © T+ and T := T~ © T+. Similarly, A^ is constant with 
respect to 8 e iff (Ao , Ai ) € ker (t _ © 5 + © 5~ © t + ) . However, it is easy to see that the existence of a nontrivial 
<5 e -constant spinor field implies the vanishing of the curvature of d e , and then, via e.g. the Gauss-Bonnet 
theorem, that the 2-surface must be a torus (see [50]). Thus if we want to weaken the condition S c Xa = 0, 
we naturally arrive at the first order operators S := S + © <5~, h~ := S~ © h + := S + © t + , c+ := S + © t~ , 
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cT := 5~ ® t+ and t := t~ ® t + . (This 8 should not be confused with the operator S := m a \7 a of the 
Newman-Penrose formalism introduced in section 3.) 

For round spheres of area-radius r (i.e. for a metric sphere in a spherically symmetric spacetimc whose 
radius r is defined by Airr 2 := Area(S)) one has a = = er', the convergences p and p' are constant on S, 
and d — \ od and d' — \ od' . In particular, in Minkowski spacetime p = — i and p' = ^. By (3.2-7) and 
(3.9), asymptotically, near the future null infinity, the operators (A.0.1.a-4.a) tend to oc T _ (Ao, Ai) := qOXq, 
OC A+(A , Ai) := od'Xo - Ai, - oo A"(A , Ai) := odXi + ±A and -ooT + (X , Ai) := 0$% - ct°A , respectively. 
The analogous oo^* and are just the corresponding unit sphere edth and edth-prime operators. In the 
rest of this appendix we discuss the direct sum operators briefly, calculate their kernels explicitly for round 
spheres, and determine how they are related to the asymptotic symmetries of the spacetimc at the future 
null infinity. Although at the quasi-local level only TC ± yield acceptable spinor fields for (2.2) and (2.3) [47], 
we will see that asymptotically, near the future null infinity any of h~ , Hr , C + , T and t (and 'at' infinity 
ocA also) can be used to recover the asymptotic symmetries. 

A.l The Dirac Witten operators A and 5 

By the definitions Aa>aX a = M' A+(A , Ai) — 5a> A~ (A , Ai) and 5a<aX a — tA"5 + (Ao,Ai) — oa>8~(\o, Ai); 
i.e. A(Ao, Ai) = and <5(Ao, Ai) = are just the Dirac- Witten equations on S with respect to A e and S e , 
respectively. They are elliptic operators with vanishing index, and hence in the generic case their kernel is 
zero dimensional. In fact, by dimker5(_! ) = dimkerc)'^ ) — one has dimker<5 = 0. On the other hand, 
there might be exceptional 2-surfaces, even among the round spheres, for which dimker A is not zero. To 
find these exceptional round 2-sphcres first take the 5-derivative of A+(A ,Ai) = and the e*'-derivative 
of A~(Ao,Ai) = 0. We obtain dd'Xo = pp'^a and d'd\\ = pp'Ai, i.e. Ao and Ai are eigenfunctions of 
dd' and d'd, respectively, with the same eigenvalue pp' . Thus let us expand them as series of the s = ±| 
spin weighted spherical harmonics: A = J2j. m Co™ an d ^1 = J2j, m ^i™ -i^},m, where j = |,§,-- 

and to = + (see e.g. [40]). Recalling that d s Yj,m = \f [j + s + - s) s +iY 3 . m and 

S>' sYj. m = -^,\f(j — s + + s) s _iYj tTn , the second order equations yield that 

2rV = -(, + ^) 2 , 3 = \\- (Al.l) 

i.e. the Dirac-Witten operator A can have non-trivial kernel only for discrete values of2r 2 pp'. Thus for a 
given, allowed pp' labelled by j one has A = J2 m c o™ i^S'.m an d ^1 = Em C T -±Yj,m- Substituting these 
into the first order equations we finally get c™ = y/2(j + ^)~ 1 rp' c™, i.e. 




(A1.2) 



In particular, for j — \ (A. 1.1) is just the condition that the round sphere is a metric sphere in Minkowski 
spacetime, whenever (A. 1.2) is the combination of the restriction to S of the two constant spinor fields given 
explicitly by (3.1). (For the explicit expression of s Yj^ m , see e.g. [40,54].) However, apart from the Minkowski 
case, (A. 1.1) yields negative Hawking energy E H := (Area(5 r )/167rG 2 )5(l + ^. § pp'AS r ) = ^(1— (i+5) 2 ), 
which is just the (holomorphic and anti-holomorphic) Dougan-Mason energy in an appropriate basis, because 
the mass is to 2 := £ab£A'B'-P AA P bb = ^y(l + 2r 2 pp') 2 . Thus, in general, the operators A and 5 do not 
define special spinor fields on round spheres. 

On the other hand, on large spheres near the future null infinity the Dirac-Witten equations tend 
to ooA(Ao,Ai) = 0, the Dirac-Witten equations on the unit round sphere in Minkowski spacetime. But 
oo A(Aq, Ai) = admits £q = A and £f = I A , given explicitly by (3.1), as independent solutions. Therefore, 
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by the analysis of Section 3, both the BMS translation and (up to supertranslations) the BMS boost-rotation 
generators can be recovered from the (exceptional) solutions of the Dirac -Wittcn equations on the cuts of 
I + . However, apart from the Minkowski case, these solutions are not limits of solutions of the Dirac-Witten 
equations on finite, large spheres, because the latter's do not exist in general. 

A. 2 The anti-holomorphy operators H~ and h~ 

By the definitions m e A e \A — — o^A _ (Ao, Ai) — t^T^Ao, Ai) and m e 5 e \A = —oaS~(Xo, Ai) — t^i^Ao, Ai); 
i.e. Xa is anti-holomorphic with respect to A e or S e iff W~(A ,Ai) = or /i~(A ,Ai) = 0, respectively. 
The symplectic scalar product of any two anti-holomorphic spinor fields (with respect to either A e or 5 e ) is 
anti-holomorphic, and hence constant on<S. liT and h~ are elliptic, and their index is 2(1 — g), where g is the 
genus of S. Therefore, in the generic case on topological 2-spheres their kernel is two complex dimensional. 

In fact, by dimker5(_ l!0 ) = and dimker5( 10 ) = 2 the <5 e -anti-holomorphic spinor fields have the 
form \a = — Aoiyt, where Ao € ker<3(i :0 ). Therefore, the space of <5 e -anti-holomorphic spinor fields does not 
inherit a natural SL(2, C) scalar product from eab- However, on round spheres these spinor fields can be 
normalized with respect to each other: For ctA, &a € kcr/i~ the combination ao<5o' + <5oao' is constant on 
S, and can be required to be y/2. In fact, on round spheres for the two independent explicit solutions we 
can choose 

n i\/2 C , i</2 
a A = —===i A , a A = = =i A , (A2.1) 

Vi + CC Vi + CC 

which arc normalized in this sense. The transformations leaving this normalization invariant is SL(2, C) x 
U(\). (In fact, the normalization condition is only one real condition, and this leaves an unspecified phase 
in the basis solutions (A. 2.1).) The spinor components a\o A are proportional to the s = | spin weighted 
spherical harmonics: They are — ^2\/2ttiYi ± i, respectively, and their normalization with respect to the 
unit sphere volume form is given by § s {a^o A ){a^,d B )d<S = 2\p2-n diag(l, 1) = 47rer AB . Since on large 
spheres the equations for the <5 e -anti-holomorphic spinor fields are just those on the round unit sphere, 
the contravariant form of the independent (5 e -anti-holomorphic spinor fields on cuts of X + are given by 
a± := e AB a%e BA ; i.e. explicitly by a A = — i</2(l + CC)"^ and a i = iv / 2C(l + (()~^ A - Comparing 
these with the functions ta at the end of the first paragraph in Section 3, one can see that aj^ — t^l a — ta.I A , 
which are just the leading parts of the spinor constituents of the BMS translations at I + . The main part of 
the anti-self-dual BMS rotations can also be expressed as — OAOB^^a^y These <5 e -anti-holomorphic spinor 
fields are used in [55] to find a relationship between the Bondi-Sachs mass at the past null infinity and the 
area of a marginally trapped surface. 

For the independent A e -anti-holomorphic spinor fields on round spheres we choose 

^ = -?17TTcc-( 2 ^-^)' ^-T27^T^ OA + vb LA ^ (A2 - 2) 

These are normalized with respect to the pointwise SL(2, C) scalar product: s AB v^v^ — e AB . Therefore, a 
natural spin space structure is inherited on the space of A e -anti-holomorphic spinor fields on round spheres. 
The L 2 scalar product of the spinor components are j> 5 ^ A ^ B 'dS = 2\/27rdiag(l, (2rp')- 2 ) and § s v^Vy dS = 
\f2n diag((2rp') 2 , 1). In Minkowski spacetime v A reduces to the restriction to S of — £ A = {I a, — Oa}, minus 
the dual of the Cartesian spin frame (3.1). Since the equations for the A e -anti-holomorphic spinor fields on 
large spheres are the ones on the round unit sphere in Minkowski spacetime, the spinor constituents of the 
BMS translations on the cuts ofl + are given by the components of the A e - anti-holomorphic spinor fields via 
ta = °a£a = °A£ AB VgeBA, and hence the spinor constituents themselves are £ Ab v%£-ba — tal a + 0(Q). 
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This representation of the BMS translations is used in [35]. The expression of the main part of the anti-self- 
dual BMS rotations by the A e -anti-holomorphic spinors is then obvious. 



A. 3 The holomorphy operators H + and h + 

As a simple consequence of the definitions, A^ is holomorphic with respect to A e or to 5 e iff H + (\o, Ai) = or 
h + (Xo, Ai) = 0, respectively. The properties of these differential operators are similar to those of H~ and h~ , 
thus we concentrate only the differences. By dimkerc5( 10 ) = and dimker5(_ 10 ) = 2 the 8 e -holomorphic 
spinor fields have the form A^ = AiOa, where Ai £ kercJ(_ 10 ). On round spheres for the two independent 
explicit solutions we can choose 

whose components (3^t A are proportional to the s = — \ spin weighted spherical harmonics: They are 
— \/2 y/lr _ i Yi ± i , respectively. Their contravariant form are the spinor constituents of the BMS translations 
on the past null infinity Z~ , and not on X + . On round spheres 



^A = 4m^=f(-°A + CV2i A ), p\ = 4m -j^^Uoa-tpVIla) (A.3.2) 
V2 y/l + (( Vrp J V2 Vl + CC V ' 

form a normalized spin frame in the space of A e -holomorphic spinor fields, which in Minkowski spacctime 
reduce to {I a, — Oa}- However, the A e -holomorphic spinor equations on large spheres near the future null 
infinity do not reduce to the ones on a round sphere: Because of the <t° term in ryo T + the solutions of 
oo 7i + (A , Ai) = are not the restriction to S of the constant spinor fields of the Minkowski spacetime. They 
are constant only if &° = 0, i.e. in absence of outgoing gravitational radiation. 

A. 4 The operators C + and c + 

By the definitions the equations C ± (Ao,Ai) = and ^(AcAi) = can be written into the manifestly 
covariant form (A a As)7r ±B c = and ((5 a As)7r ±s ' c = 0, respectively, where tt +A b = l A ob and -k~ A b = 
— o A lb are the projections of the spin spaces onto the space of the ±1 eigenspinors of the spinor "f A b = 
o a ib + t A os, defining a GHP-spin frame independent notion of chirality on the spin spaces (see [51]). C ± 
and c ± are not elliptic operators. In fact, by dim ker 5(^0) = 2 and dimkercJ^ ) = on round spheres 
the kernel of C + is two dimensional only if p ^ 0, otherwise dimkerC + = oo. In particular, c + is just the 
p = special case of C + , whenever the elements of c + have the form A^ = \\oa with arbitrary Ai : S — > R. 
On round spheres with non-zero p the kernel of C + is just the kernel of TL + ; i.e. for the two explicit 
solutions we can choose (A. 3. 2). However, the equations C + (Ao,Ai) = on large spheres near the future 
null infinity tend to the asymptotic twistor equations of Bramson [53], which turn out to be the equations 
C + (A , Ai) = on the round unit sphere in Minkowski spacetime, and hence kerooC + = kcr oc W + still holds. 
Therefore, the l a -components of the normalized solutions of the asymptotic twistor equation reproduce the 
components ta of the spinor constituents of the BMS translations on I + , and the solutions themselves are 
s AB p^Eba = tal a + 0(fl). This is one of the most popular representation of the BMS translations (see 
[22-24,34,38,40,56]). The expression of the BMS rotations is then straightforward. 

A. 5 The operators C and c~ 

By dimkerc^! o) = and dimkcrc?(_ 1 ) = 2 the kernel of C~ on round spheres is two dimensional only if 
p' 7^ 0, otherwise dimkerC" = oo. c~ is the p' = special case of C", and the elements of c~ have the form 
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\a = — A /-a with arbitrary A : S — > R. If p' ^ then kerC~ = kcrW~, thus for the two explicit solutions 
of C~(Ao, Ai) = we can choose (A. 2. 2). However, the equations C~(\q, Ai) = on large spheres near the 
future null infinity are not the equations C~(\q, Ai) = on the round unit sphere in Minkowski spacetime, 
just because of the presence of the a term in oc^" + - Therefore, the solutions of C~(A , Ai) = on large 
spheres can be used to represent the BMS translations on the past null infinity. 

A. 6 The 2-surface twistor operators T and t 

The 2-surface twistor equations T(Ao, Ai) = can be written into the covariant form Te'ea B ^b '■= ^-e'(e^a) 
+ \lEAl BC Ae>b^c = 0; i.e. the vanishing of the 7 J 4B-trace-free symmetrized A e -derivative of \a (see [45]). 
Similarly, i(Ao,Ai) = is equivalent to 5e<(e^a) + \lEAl BC '°~e' b^c = 0. T and t are elliptic operators 
with index 4(1 — g). On round spheres they coincide, and by dimkerc^o) = dimkere^^ ^ = 2 their kernel 
is four dimensional. Therefore, the general explicit solution of the 2-surface twistor equations has the form 
ua = sa u ^ + oaPa f° r complex constants <za, &a, where a A and (3 A are given by (A. 2.1) and (A. 3.1), 
respectively. Introducing ita> '■= i Aa>a^ A and the 2-surface twistor Z a :— (u) A ,tta'), the usual Hermitian 
pointwise scalar product of any two 2-surface twistors Z" and Z ,a is constant on S, and hence defines the 
familiar Hermitian scalar product on the space T Q of the 2-surface twistors. However, the flat spacetime 
definition of the infinity twistor docs not define any (global) twistor on T", because e A B tta'Tt' B / is not 
constant on S. Its modification, I a/ 3Z a Z'^ :~ e A B 'it a'^'b' + 2^(1 % r2 PP') £ abu a lo' b is constant on S, but, 
apart from the Minkowski case 2r 2 pp' = — 1, its rank is four and hence it fails to be simple. Thus, even 
on general round spheres, Iq^Z' 3 = cannot be used to define 'translations' in the twistor space T Q . The 
primary parts of the 2-surface twistors defined by t and T coincide, but their secondary parts do not. The 
2-surface twistors defined by t can be recovered from those of T as the p = = p' special case. 

On large spheres near X + the solutions of the 2-surface twistor equations have the structure loa = 
aA&A + °aPa + a Al A °A = (°c £CBa B a A-) L A + {°aPb lB + a Al A )°Ai where 7 A are of spin weight — ^ and 
satisfy o$'7 A = 5-°a A . Comparing this with the results of subsection A. 2 we find that lu a = ap>_T A L A + 
(&A/3 A + aA7 A )o A = aAT A L A + 0(Q); i.e. the primary part of the 2-surface twistors Z a — (u> A , i Aa'b^ B ) are 
just the spinor constituents of the BMS translations. Although the 2-surface twistors on large spheres still 
form a 4 dimensional complex vector space, half of them (parametrized by the 6a 's) die off asymptotically. 
In absence of outgoing gravitational radiation ooT reduces to the homogeneous ooi, whose solutions are just 
those of T and of t on the unit round sphere. Interestingly enough, in the conformal approach of the 2-surface 
twistor equation on the cuts of I + (see [12-20,38,40]) the absence of outgoing gravitational radiation yields 
the homogeneous equations o$A = 0, od'\\ = only after an appropriate supertranslation. 
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